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Abstract 

We prove a general theorem on spectral convexity with respect to particle number for 2N de- 
generate components of fermions. The number of spatial dimensions is arbitrary, and the system 
may be uniform or constrained by an external potential. We assume only that the interactions are 
governed by an S'Lf (2iV)-invariant two-body potential whose Fourier transform is negative definite. 
The convexity result implies that the ground state is in a 2iV-particle clustering phase. We discuss 
implications for light nuclei as well as asymmetric nuclear matter in neutron stars. 
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Interacting fermions with more than two components exhibit a variety of low tempera- 
ture phenomena. Of particular interest are phenomena which appear in different quantum 
systems and therefore could be characterized as universal. One example in three dimen- 
sions is the Efimov effect, which predicts a geometric sequence of trimer bound states for 
interactions in the limit of zero range and infinite scattering length The 
Efimov effect is forbidden for two-component fermions due to the Pauli exclusion principle 
but can occur for more than two components. Efimov trimers have recently been observed 
in ultracold cesium as indicated by a large three-body recombination loss near a Feshbach 

Q 

resonance [7[ . Once the binding energy of the trimer system is fixed for interactions at zero 
range and large scattering length, the binding energy of the four-body system is also deter- 
mined. This is in direct analogy with the Tjon line relating the nuclear binding energies of 
3 H and 4 He Q, [l0|. In two dimensions a different geometric sequence has been predicted 
for zero-range attractive interactions. In this case the geometric sequence describes the 
binding energy of iV-body clusters as a function of N in the large N limit U, 12, Q]- 



Several recent studies have investigated pairing and the superfluid properties of three- 



component fermions 15j, [16|, [17|, [18]. Systems involving four-component fermions are of 
direct relevance to the low-energy effective theory of protons and neutrons. Due to an- 
tisymmetry there are only two S-wave nucleon scattering lengths, corresponding with the 
spin-singlet and spin-triplet channels. Some general properties of this low-energy effective 



;heory have been studied such as pairing, the 
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22j . Wu and collaborators 
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ermion sign problem, and spectral inequalities 



251 ] have pointed out that the effective theory 



has an accidental 50(5) or Sp(A) symmetry, and several different phases such as quintet 
Cooper pairing or four-fermion quartetting could be experimentally realized for different 



scattering lengths with ultracold atoms in optical traps or lattices [26|, [271]. When the 
scattering lengths are equal the symmetry is expanded to SU(4). This symmetry was first 
studied by Wigner and arises naturally in the limit of large number of colors for quan- 
tum chromodynamics j^. fsol] . The fact that both the spin-singlet and spin-triplet nucleon 
scattering lengths are unusually large means that the physics of low-energy nucleons is close 



to the Wigner limit 3l|, | 

In the following we prove a general theorem on spectral convexity with respect to particle 
number for 2N degenerate components of fermions. The theorem holds for any number 
of spatial dimensions, and the system may be either uniform or constrained by an external 
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weak attractive potential 
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FIG. 1: Illustration of the convexity constraints for the ground state energy E as a function of 
particle number A. The line segments show the convexity lower bounds. 

potential. We assume only that the interactions are governed by an S'?7(2iV)-invariant two- 
body potential whose Fourier transform is negative definite. The main result is that if the 
ground state energy E is plotted as a function of the number of particles A, then the function 
E(A) is convex for even A modulo 2N. Furthermore E(A) for odd A is bounded below by 
the average of the two neighboring even values, E(A — 1) and E(A + 1). This is illustrated 
in Fig. [T] for both the weak attractive and strong attractive cases. This convexity pattern 
could be regarded as an SU(2N) generalization of even-odd staggering for the ground state 
energy in the attractive two-component system. 

A weaker form of this inequality was proven for A < 2N and zero-range attractive 
interactions [21]. Here we extend the proof to any A and any S , f/(2A^)-invariant potential 
with a negative-definite Fourier transform. Another difference between this and the previous 
analysis is that we use a fixed particle number formalism rather than the grand canonical 
formalism. This is essential for strongly-attractive interactions and the requirement of 
keeping the physics in the low-energy regime. By limiting the number of particles we avoid 
a collapse towards high densities and the need for hard core repulsion for stability. 

We start by considering 2N degenerate components of nonrelativistic fermions in d spatial 
dimensions. We assume the interactions are governed by an S77(2AQ-invariant two-body 
potential V{f) whose Fourier transform V(p) is strictly negative. We also allow an SU(2N)- 



invariant external potential U{r) whose properties are not restricted. The general form of 
the Hamiltonian is 

# = --*- I ' d d fa\{r)V 2 a l {r)+ I ' d d fU(f)p(r) 

m i=l,-,2N •* •* 

+ i J d d rd d ? : p(r)V(f- f)p(f) : , (1) 
where p(r) is the SU (2iV)-invariant density, 

p(o= E ^xo- (2) 

i=l,- ,27V 

The : symbols denote normal ordering. We consider the system on a hypercubic lattice 
using a transfer matrix formalism. We let n = (n s ,n t ) represent (d+ l)-dimensional lattice 
vectors. The subscript s on n s denotes a d- dimensional spatial lattice vector. We write 
the c?-dimensional spatial lattice unit vectors as 1, • • • , d. Throughout our discussion of 
the lattice system we use dimensionless parameters and operators which correspond with 
physical values multiplied by the appropriate power of the spatial lattice spacing a. We let 
a t be the temporal lattice spacing and a t be the ratio a t /a. L denotes the spatial length of 
the periodic hypercubic lattice. 

We use the notation V(2nk s /L) for the Fourier transform of the lattice potential V(n s ), 

V{2nk s /L) = V{n s )e i2 ™ a ° /L . (3) 

n s 

~ —* 

By assumption V(2irk s /L) is strictly negative. Let M be the normal-ordered transfer matrix 
operator 



M = : exp 



-a t H hee - a t ^2,U{n s )p{n s ) - y ^ p(n s )V(n s - n' s )p{n' s ) 



(4) 



where H^ ee is the free lattice Hamiltonian, 

Hfrcc = ~2m ^ E { a l(^ s ) ai(n s + L) + ai(n s - l s ) - 2ai(n s )^ . (5) 

n s L=i,-,di=i,-,2Af 

Let V~ l {n s ) be the inverse of V(n s ), 

i p —i2nn a -k s /L 

L d f V(2nk s /L) 



We now rewrite powers of M using an auxiliary field 0, 

M Lt = jD(t)e~ sw M L( _i(0) x ••• x M o (0), (7) 

where 

S (<f>) = "f E E -<)0«,n t ), (8) 



2 

ret fts ,n' s 



M nt (4>) = : exp 



-a t H fTCC -a t y^U (n s )p{n s ) + a t E Wt)p(«« 



(9) 



fe s n 3 ,n t V / * 

Let f^(n s ), f^(n s ), ■ ■ ■ be a complete set of orthonormal real- valued functions of the 
spatial lattice sites n s . We refer to these functions as orbitals. We take f^(n s ) to be 
strictly positive but otherwise regard the form for the orbitals to be arbitrary. If the total 
number of lattice sites is L d then we have a total of L d orbitals. We denote a one-particle 
state with component i in the k th orbital as f- k ^. 

Let B and C be any finite subsets of the orbital indices, B, C C {1, 2, • • • , L d }. From these 
we define \&C 2N ~^) as the quantum state where each of j components fill orbitals B and 
each of the remaining 2N — j components fill the orbitals C. The order of the component 
labels is irrelevant, and so we assume that the first j components fill orbitals B and last 
2N — j components fill orbitals C. The total number of fermions in state \B j C 2N ~ j ) is 
j \B\ + (2iV — j) \C\ , where \B\ and \C\ are the number of elements in B and C respectively. 

We define E&j&N-j as the energy of the lowest energy eigenstate with nonzero inner 
product with \B j C 2N ~ j ). We let Z^) Q2N _ 3 be the expectation value of M Lt for \B j C 2N ~ j ), 

Z^ C2N _, = (B j C 2N ~ j \ M Lt \B j C 2N ~ j ) . (11) 



In the limit of large L t the contribution from the lowest energy eigenstate dominates and 
therefore 

lv, ( 7 L t ^ 

We can write Z^ t JC2N _ J using the auxiliary field 0, 



E &C 2N- 3 = - lim !^g£^ziZ . (12) 



Z L B y N ^ = jD<p e~ s ^ (B j C 2N ~ j \ M Lt _x(0) x • • • x M o (0) \B j C 2N ~ j ) . (13) 
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At this point we define matrix elements for the one particle states, 



Mvjktt) = (fP\ M^M) x ■ ■ ■ x M o (0) \f™ 



(14) 



The component index i in (fl4"|) does not matter due to the SU(2N) symmetry. Each entry 
of the matrix M.k',k{4>) is real. We let Mb{4>) be the \B\ x \B\ submatrix consisting of the 
rows and columns in B and let Aic{4>) be the \C\ x \C\ submatrix for C. 

Each normal-ordered transfer matrix operator M nt (0) has only single-particle interactions 
with the auxiliary field and no direct interactions between particles. Therefore it follows 
that 



-S(4>) 



[detM B (<j>)) j [det M c {(j))] 



2N-j 



(15) 



This result is perhaps more transparent if we pretend for the moment that each of the j \B\ + 
(2N — j) \C\ particles carries an extra quantum number which makes them distinguishable. 
We label the extra quantum number as X. So long as the initial and final state wavefunctions 
are completely antisymmetric in X for particles of the same component then this error in 
quantum statistics has no effect on the final amplitude. So we can factorize the transfer 
matrices M nt (<p) as a product of transfer matrices for each X. This then leads directly to 

(US}. 

Let ni an d be integers such that < 2ni < j < 2ri2 < 2N. Let us define the new 
positive-definite measure 



Dcj) = D<f> e~ sw [&et M B ((f))Y ni [detM c (4>)} 



2ni 



i2jV-2n 2 



(16) 



so that 



J B j C 2N-j 



D<t> [det Mb 



Then by the Holder inequality \Z 



j-2m 



[det M c (<f>)\ 



2n2— j 



BiC 2N -i 



is bounded above by 



D(j) \detM B ((j>)\ ' 



2ri2— 2ni 



D(f) \det M c (<i>)\ 



2n 2 — 2ni 



( r/Lt \ 2n 2 -2n 1 ( N 2n 2 -2n x 

jg2n 2 c2N-2ri2 ) \^ g2n 1 Q2N-2n 1 j 



(17) 



(18) 



Taking the limit L t — > oo we deduce that the energies satisfy the inequality 



Er- 



I \, j — 2ni 7-1 ZU2-] 771 

BJC^-j > 2n2 _ 2ni -£j B 2n 2C 2N-2n 2 + 2n 2 — 2ni -"-B 2 ™lC 2JV_2rl l 



2n2 — j 



(19) 
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This is a statement of convexity for E B j C 2N-j as a function of j between even endpoints 
j = 2n\ and j = 2n2- If we now take \B\ — K + 1 and \C\ = K, then the total particle 
number is A = 2NK + j and A lies between 2NK and 2N(K + 1). The inequality in (TT9"j) 
is precisely the convexity pattern in Fig. [TJfor -£■(^4) as a function of particle number. 

We point out that for the special case K = 0, we can take B to be the first orbital 
and C to be the empty set. In this case Aisi^) is simply a number. Furthermore since 
f^(n s ) is strictly positive, M.b{4>) is also positive so long as the temporal lattice step a t 
is not excessively large. Since det M.ei.'P) — -Mb(0) > it is no longer necessary that the 
power of det -Mg(0) be even to insure positivity. Therefore E(A) is actually convex for all 
A between and 2N and not just even A. 

These convexity relations could be checked using any number of attractive SU(2N) mod- 
els in various dimensions. This will be checked in future studies. Here we instead examine 
actual nuclear physics data to investigate Wigner's approximate SU{4) symmetry in light 
nuclei. It is by no means clear that the interactions of nucleons in light nuclei can be 
approximately described by an attractive 5*?7(4)-svmmetric potential. Recent results from 
nuclear lattice simulations hint that this might be possible 33j, |34|, however there are forces 
even at lowest order in chiral effective field theory which break SU(4) invariance in addition 
to being repulsive. Nevertheless all of the SU(4) convexity constraints are in fact satisfied 
for the most stable light nuclei with up to 16 nucleons as can be seen in Fig. [21 The line 
segments drawn show all of the convexity lower bounds. 
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here have been several recent studies of alpha clustering in nuclear matter 



39 



401 ] as well as multiparticle clustering in other systems 24, 



2a 



26 



27|, 



41 



35 



42 



36 



37|, 



The 



results presented here give sufficient conditions for the onset of this multiparticle clustering 
phase. One can also make a definite prediction about the j-component quasiparticle energy 
gaps. Starting from a 2iVi^-fermion S'L r (2A^)-symmetric state, let 5j be the extra energy 
required per fermion to add j fermions, all of different components. The ground state energy 
for 2NK + j fermions is a convex function for even j in the interval from j = to j = 2N. 
Therefore it follows that 82 > ^4 > • • • > 52N- Since the ground state energy for 2NK + j 
fermions is also convex for j = 0,1,2, we conclude furthermore that 5i > 8 2 > 64 ■ ■ ■ > S 2 n- 
We note that for the strongly attractive case these energy gaps are negative, and it is more 
natural to speak of energy gaps per missing fermion for the corresponding j-component 
quasiholes, 6j. In this case we find again > 5% > 5% ■ ■ ■ > 5% N . 
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FIG. 2: Plot of the energy versus particle number for the most stable light nuclei with up to 16 
nucleons. The line segments show the convexity lower bounds. 

In summary we have derived a general result on spectral convexity with respect to particle 
number for 2N degenerate components of fermions. We assume only that the interactions 
are governed by an SU (2iV)-invariant two-body potential whose Fourier transform is neg- 
ative definite. The ground state energy £ as a function of the number of particles A is 
convex for even A modulo 2N. Also E(A) for odd A is bounded below by the average of 
the two neighboring even values, E(A — 1) and E(A + 1). When applied to light nuclei 
for A < 16 all of the convexity bounds for 577(4) are satisfied. These results give further 
evidence that an approximate description of light nuclei may be possible using an attractive 
SU (4)-symmetric potential. This would be a direction worth pursuing since the same theory 
could then be applied to dilute neutron-rich matter with a finite number of protons. The 
residual 577(2) x £77(2) symmetry for proton spins and neutron spins would guarantee that 
the Monte Carlo simulation could be done without fermion sign oscillations. The physics of 
this quantum system would be helpful in understanding the superfluid properties of dilute 
neutron-rich matter in the inner crust of neutron stars. 

The author thanks Gautam Rupak and Thomas Schdfer for discussions and correspon- 
dence with Philippe Lecheminant and Congjun Wu. This work was supported in part by 
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